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r"| ■ NLO-QCD corrections to vector boson pair production via weak boson 

f> { fusion have recently been calculated and implemented into flexible parton- 

level Monte-Carlo programs. These allow for the computation of cross sec- 
tions and kinematical distributions within realistic experimental cuts. We 
r~| ', summarize the basic elements of the calculation and review phenomenolog- 

ical results for the LHC. 
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1. Introduction 

One of the major goals of the CERN Large Hadron Collider (LHC) is 
the understanding of the mechanism of electroweak (EW) symmetry break- 
ing pQ. In this context, vector-boson fusion (VBF) processes form a promis- 
ing class of reactions: Higgs production in VBF, i.e. the reaction qq — > qqH , 
has been proposed as a particularly clean channel for the discovery of the 
Higgs boson and a later determination of its couplings [2J. An important 
background to the H — * VV decay channel (V = W or Z) in VBF is 
constituted by continuum VV production via VBF, i.e. EW pp — ► VVjj 
production The precise knowledge of the standard model cross section 
for this reaction at NLO-QCD accuracy becomes crucial for distinguishing 
enhancements in VBF processes due to signatures of new physics, such as 
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Fig. 1. Representative diagrams contributing to the LO cross section for the process 
qq — > qq£ + £~f + £'~ (a) and to the virtual corrections for the same reaction (b). 

strong electroweak symmetry breaking from possible effects of higher 
order perturbative corrections. 

We have therefore calculated the NLO-QCD corrections to the reactions 
pp -> e+Ve^v^jj 5 ], pp -> e+e" fi+fi'jj, and -> e + e~ v^v^jj 0. In 
the following we will refer to these processes briefly as "EW VVjj pro- 
duction" . Our calculations have been turned into fully flexible parton-level 
Monte-Carlo programs allowing for the computation of cross sections and 
distributions within realistic experimental cuts, in complete analogy to the 
similar programs for the Hjj signal [7] and Vjj production 8 via VBF. 
In this contribution we briefly recollect the elements of our calculation and 
discuss the basic features of our results. 

2. Elements of the Calculation 

The computation of partonic matrix elements for EW VVjj production 
is based on the amplitude techniques of Ref. Special emphasis is put 
on the development of a numerically stable and fast code. That is achieved 
by organizing the calculation of the matrix elements in a modular way such 
that building blocks that are encountered in several diagrams are evaluated 
only once per phase space point, as described in some detail in Refs. 011)]. 
This feature becomes particularly important for the computation of the 
real emission corrections to VVjj production due to the large number of 
contributing diagrams. 

Singularities emerging from soft and collinear configurations are regu- 
larized in the dimensional-reduction scheme ^U] with space-time dimension 
d = 4— 2e. The cancellation of the divergences with the respective poles from 
the virtual contributions is performed by introducing the counter-terms of 
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the dipole subtraction method The analytical form of the phase-space 
integrated subtraction terms, obtained after the factorization of the parton 
distribution functions, is given by 
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in the notation of Ref . , with Ai b denoting the amplitude of the corre- 
sponding Born process and Q 2 the momentum transfer between the initial 
and final state quark in Fig. ^ The computation of the virtual corrections 
requires the evaluation of self-energy, triangle-, box-, and pentagon contri- 
butions on either the upper or the lower quark line, as sketched in Fig. Q (b) 
for one representative diagram. Contributions from graphs with gluons at- 
tached to both the upper and lower quark lines vanish at order a s . Putting 
all virtual contributions together, we find 
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with Cvirt = tt 2 /3 — 7, and a completely finite remainder My. The diver- 
gent pieces in this expression exactly cancel the poles of the counter-terms 
in Eq. (|2.1jl . The remaining integrals are finite and can be computed nu- 
merically in d = 4 dimensions by means of a Passarino-Veltman tensor 
reduction ^2] ■ Numerical stability is achieved by the repeated use of Ward 
identities which allow us to express a large fraction of the pentagon contri- 
butions by a combination of box-type diagrams. 



3. Predictions for the LHC 

The cross-section contributions discussed above are implemented in fully 
flexible parton level Monte-Carlo programs which allow us to calculate cross 
sections and kinematical distributions for EW WWjj and ZZjj production 
at NLO-QCD accuracy within typical experimental acceptance cuts. We use 
the CTEQ6M parton distributions with a s {mz) = 0.118 at NLO, and the 
CTEQ6L1 set at LO 13 . We chose m z = 91.188 GeV, m w = 80.419 GeV 
and Gf = 1.166 x 10 _5 /GeV 2 as electroweak input parameters, from which 
we obtain o>qed = 1/132.54 and sin 2 9w = 0.22217. Jets are reconstructed 
from final-state partons employing the kx algorithm j!41 115j with resolution 
parameter D = 0.8. Throughout, we set fermion masses to zero and neglect 
external b- and i-quark contributions. 
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In the following, we consider EW ZZjj production within generic cuts 
that are relevant for VBF studies at the LHC. We require at least two hard 
jets with 

p Tj >20GeV, \y 3 \ < 4.5, (3.1) 

where yj is the rapidity of the (massive) jet momentum which is recon- 
structed as the four- vector sum of massless partons of pseudorapidity \rj\ < 
5. The two reconstructed jets of highest transverse momentum are called 
"tagging jets" . To suppress backgrounds to VBF we impose a large rapidity 
separation of the two tagging jets, 

&Vjj = \Vji - Vh\ > * • (3.2) 

In addition, we require the two tagging jets to reside in opposite detector 
hemispheres, 

Vh x Vh < » ( 3 - 3 ) 

with an invariant mass 

M j:j > 600 GeV , (3.4) 

and adopt the lepton cuts 

pre > 20 GeV , \r] £ \< 2.5 , AR j£ > 0.4 , 

m u > 15 GeV , AR U > 0.2 , (3.5) 

where ARj£ and Ai2^ denote the jet-lepton and lepton-lepton separation 
in the rapidity-azimuthal angle plane, respectively, and ma the invariant 
mass of an e + e~ or n~ pair. 

The total cross section for EW VVjj production contains contributions 
from the Higgs resonance with H — > VV decays, as well as from VV con- 
tinuum production. In the following we will focus on the VV continuum by 
imposing a cut on the four-lepton invariant mass 

M V v = y/ipt! +Pe 2 +P£ 3 +Pi 4 ) 2 > m H + 10 GeV , (3.6) 

where the denote the four-momenta of the leptons produced in the spe- 
cific reaction under consideration. 

The total continuum cross section for the reaction pp — > e + e~ at 
NLO within the cuts of Eqs. (|3.1H3.6j) and a Higgs mass of mu = 120 GeV is 
displayed in Fig. 21 as a function of the renormalization scale [Ir = ^mz with 
a fixed factorization scale \ip = mz (dashed green curve), as a function of 
Hp = fymz with fj,R = mz (dot-dashed blue curve), and for the case, where 
both scales are varied simultaneously, fin = hp = £,mz (solid red curve). 
The LO cross section depends only on fip = £,mz (dotted black curve). 
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Fig. 2. Dependence of the total pp — > e + e jj continuum cross section at 

the LHC on the factorization and renormalization scales. 



As in the analogous cases of e + v e ^i~9^ and e + e~ v^v^ production via VBF, 
discussed in Refs. jSHH], the scale variations of the NLO prediction are found 
to be below the 2% level when the scale parameter runs from £ = 0.5 to 
£ = 2.0, while the LO cross section changes by more than 20% in the same 
range of £. 

Figure |31 shows the transverse momentum distribution of the tagging jet 
with the highest pr for EW e + e~ fi + \i~ jj production. The factorization 
and renormalization scales are fixed to [if = /j>r = mz- The K factor, 
defined by 

= da NLO /d X 

and shown for this distribution in Fig. 01(b), illustrates the change in shape 
when going from LO to NLO. At NLO, smaller values of are preferred by 
the tagging jet, which is due to the extra parton which emerges in the real 
emission contributions. In a similar way, other distributions such as the in- 
variant mass distribution of the tagging-jet pair, da/dMjf (see Ref. 0), or 
the transverse momentum distributions of the detected leptons are affected 
by the inclusion of NLO corrections. On the other hand, some angular 
distributions such as da/dr]™ ax , where r/™ a:r designates the largest lepton 
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Fig. 3. Transverse momentum distribution of the tagging jet with the highest pt in 
EW e + e~ fi + [J.~jj production at the LHC. Panel (a) shows the NLO result (solid 
red line) and the LO prediction (dashed black line) . Panel (b) displays the K factor 
defined in Eq. (|3~7|) . 

rapidity emerging in the scattering, barely change at NLO as illustrated 
in Fig. 0J The step in da/drfl 1 ^ at r]'^ iax = 2.5 is due to the rapidity cut 
imposed on the lepton rapidities in Ea. H3.5|) . 

Our discussion on the effect which higher order QCD corrections have on 
dynamical K factors has been based on EW e + e~ production. How- 

ever, the shapes of kinematical distributions in the other VBF production 
processes we have studied, i.e., pp — > e + e~ v^v^jj and pp — > e + v e p~v )1 jj, 
behave very similarly. 



4. Conclusions 

In this contribution we have reviewed some results on EW WWjj and 
ZZjj production at NLO-QCD accuracy obtained with fully flexible Monte- 
Carlo programs. The higher order corrections to these reactions turn out 
to be under excellent control, as indicated by the small scale dependence of 
the total cross sections and K factors close to one. We found, however, that 
some kinematical distributions exhibit a noticeable change in shape when 
NLO corrections are considered. This indicates the importance of including 
NLO-QCD contributions in precision studies of VBF processes at the LHC. 
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